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Abstract 



A large part of the mSUGRA parameter space satisfying the WMAP constraint on the 
dark matter relic density corresponds to a higgsino LSP of mass ~ 1 TeV. We find a promising 
signal for this LSP at CLIC, particularly with polarized electron and positron beams. One 
also expects a viable monochromatic 7-ray signal from its pair annihilation at the galactic 
center at least for cuspy DM halo profiles. All these results hold equally for the higgsino 
LSP of other SUSY models like the non-universal scalar or gaugino mass models and the 
so-called inverted hierarchy and more minimal supersymmetry models. 
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1. Introduction 



Weak scale supersymmetry (SUSY) is the most popular extension of the standard model 
(SM) because it is endowed with three unique features |T]. It provides (1) a natural solution 
to the hierarchy problem of the SM, (2) a plausible candidate for the cold dark matter of 
the universe in the form of the lightest superparticle (LSP), and (3) unification of the SM 
gauge couplings at the GUT scale. However it also suffers from two problems. 

(i) Little Hierarchy Problem: The LEP limit on the mass of an SM-like Higgs boson 



requires the average top squark mass to be well above Mz [SI- In models where supersym- 
metry breaking is transmitted to the visible sector at an energy scale exponentially larger 
than the weak scale, this implies some fine-tuning of SUSY parameters to obtain the correct 
value of Mz- 

(a) Flavor and CP Violation Problem: Generic SUSY models make fairly large contributions 
to CP violating processes with or without flavor violation, as represented by the K decay 
observable ex and the fermion electric dipole moments (EDM) respectively. Predictions for 
rates of CP-conserving flavor changing processes, like fi ^ ej decays, also often exceed 
experimental limits. SUSY parameters have to be chosen carefully to control these contri- 
butions. It should be noted here that the recently advocated spht SUSY model |3j tries to 
solve the second problem by pushing up the scalar superparticle masses, but at the cost of 
dramatically aggravating the first problem. 

The minimal supergravity model (mSUGRA) is by far the simplest potentially realistic 
model of weak scale supersymmetry [S]. It provides a very economical parametrisation of 
superparticle masses and couplings on the one hand and a natural explanation of the elec- 
troweak symmetry breaking (EWSB) phenomenon on the other [6 . The model is completely 
specified in terms of four continuous parameters and one sign, namely 



The first three entries represent the universal SUSY breaking scalar and gaugino masses and 
trilinear coupling at the GUT scale, /i is the supersymmetric Higgs(ino) mass parameter, 
while tan /3 is the ratio of the two Higgs vacuum expectation values. The flavor universality of 
scalar soft breaking terms avoids problems with flavor changing processes, while CP violation 
will be under control if the parameters in Q are real. 

In going down from the GUT scale to the weak scale the SU3 x SU2 x Ui gaugino masses 
evolve like their respective gauge couplings, i.e. 



m;, > 114 GeV, 



(1) 



mQ,mi/2, Ao,tanP and sign(yu). 



(2) 



M2 



Ml 



M3 



1 



(ai/aG)mi/2 ~ (25/60)mi/^ 
{a2/aG)mi/2 ~ (25/30)mi/^ 
{a3/aG)mi/2 ~ (25/9)mi/2. 



(3) 
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The higgsino masses are simply given by the parameter. The scalar masses at the weak 
scale are also related by Renormalization Group Equations (RGE) to the GUT scale mass 
parameters of eq.©. A very important scalar mass at the weak scale is the mass of the 
Higgs boson H2 that couples to the top quark. This mass appears in the EWSB condition 



,, ., , rn?TT — rn?TT tan^ B ., 



The last equality holds for tan/3 > 5, which is favored by the LEP limit of eq.((H) |21IS1- '>^'h2 
is related to the GUT scale mass parameters by the solution of its RGE j7], 



H2 



CiiTil - C2ml/2, (5) 



where we have dropped contributions oc Aq for simplicity since they do not play any impor- 
tant role here. The coefficients 6*1,6*2 depend on the gauge and Yukawa couplings. Thanks 
to the large negative contribution from the top Yukawa coupling, C*2 — 2. On the other 
hand, |C*i| ^ 1, its value and sign depending on the exact values of SM parameters (in 
particular, on rrit and a^), on the scale msvsY where the RG evolution is terminated, and 
on tan/3, with smaller tan/3 favoring smaller (possibly negative) C*i. This makes it easy to 
obtain a negative value of mjj^ as required by the EWSB condition This is the so-called 
radiative EWSB mechanism jB]. 

Combining eqs. (jlj) and ©, one sees that /x^ is related to ml and my2 by an ellipsoidal 
equation if C*i < 0, in particular for low values of tan/3 (< 5). However, for moderate 
to large values of tan/3 (> 5), favored by the LEP limit of eq.(IT)), and large msusY, C*i 
becomes positive, leading to a hyperbolic equation. These two cases have been described 
as ellipsoidal and hyperbolic branches of mSUGRA One sees from eqs. fl3l4l5p that 
in the first case Mi < \fi\; and the lightest neutralino (LSP) is the Bino (B). However 
in the phenomenologically favored case of moderate to large tan/3(> 5) one can have both 
Ml < \fi\ and Mi > corresponding to a Bino and a higgsino LSP respectively [8jj. We 
shall concentrate in this case. Since the sign of fi is not important for our analysis, we shall 
choose only positive sign for simplicity. 

In the next section we study the dark matter (DM) relic density in the Bino and higgsino 
LSP domains of mSUGRA model following the second paper of ref.0. The relic density 
has been computed using the MicrOMEGAs code [H|. A comparison with cosmological data 
[TU] on the relic DM density shows that a large fraction of mSUGRA points satisfying this 
data come from the higgsino LSP domain with ~ 1 TeV. In the following two sections we 
study the prospects of detecting the higgsino LSP in collider and DM search experiments, 
respectively. We find good prospect of detecting this particle at a 3 TeV linear collider like 
CLIC. There is also a good prospect of detecting it in the form of TeV scale gamma ray 
line from DM pair-annihilation in the galactic center for favorable profiles of galactic DM 
distribution. In the next section we shall show that all our results hold not only in mSUGRA 
but in a host of other SUSY models as well, which can naturally accommodate a higgsino 
LSP. Finally we shall conclude with a summary of our results. 
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Figure 1: mSUGRA parameter space for rrit = 178 GeV, Aq = and tan/3 = 10 (left) and 30 
(right). The patterned regions marked I and II are disallowed by the EWSB condition and 
the constraint of a neutral LSP, respectively. In the large yellow region the LSP is Bino-like. 
Red points satisfy the constraint (8) on the dark matter relic density. 

2. Higgsino LSP as DM in mSUGRA 

Figs. 1 show the mSUGRA parameter space satisfying the EWSB condition for a moderate 
and a large value of tan f3. We have used our own code for the solution of the relevant RGE 
and the treatment of EWSB, including dominant loop corrections. The upper edge of the 
allowed (white) region corresponds to the hyperbolic boundary from eqs. ()4l5|) for the LEP 
limit of > 100 GeV 0. In fact, the bulk of this disallowed region (I) corresponds to 
/i^ < 0, i.e. no EWSB. In the bottom strip (II) the tau slepton fi becomes the LSP. This 
is disallowed by the astrophysical constraint requiring a neutral LSP f2j. Over the allowed 
region the LSP is the lightest neutralino, which can be a combination of gaugino and higgsino 
states, i.e. 

X = Xi = NuB + N,2m + iVi3^? + Ni^Hl (6) 
The gaugino component of the LSP is defined by the fraction 

Z, = N',, + N',,. (7) 

The yellow (light shaded) region in Fig. 1 corresponds to a dominantly gaugino (in fact Bino 
B) LSR 

One sees from Fig. 1 that in most of the mSUGRA parameter space the LSP is dominantly 
B. Notice, however, that the bulk of this region is disallowed by the constraint on the DM 
relic density from cosmological data, in particular from the WMAP satellite experiment, JU] 

= 0.113 ±0.017, (8) 
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where h = 0.71 ± 004 is the Hubble constant in units of 100 km/(sxMpc) and Q is the 
rehc density in units of the critical density. In fact, one usually finds an overabundance 
of DM relic density {Q^ > 1). This is because B does not carry any gauge charge and 

hence does not couple to the gauge bosons. A Bino-like LSP therefore mostly annihilates 

f - 

via the exchange of sfermions, XX — ^//^ which is suppressed by the large sfermion mass. 
Only in some special cases like ~ rrif^ and ~ mAfi can one get large co-annihilation 

X^i~^Tl and pair annihilation xx~~^bb, tt rates respectively ^]. Correspondingly one can 
see a few (red) points of acceptable B DM density near the lower boundary (co-annihilation); 
in mSUGRA resonant A— exchange becomes possible only for tan/3 > 50.^ 

Most of the points satisfying the DM relic density constraint (jH)) are seen to lie very near 
the hyperbolic boundary [9j . The few points near the lower end of this boundary correspond 
to the so-called focus point region j^, where the LSP has a significant higgsino component, 
although it may be still dominated by B. Such an LSP couples via its higgsino component 
to W and Z bosons, and through gaugino-higgsino mixing to Higgs bosons, and can thus 
annihilate into both fermionic and bosonic final states. LHC signatures for the focus point 
region have been investigated in jl4j. Note however that the large majority of DM-allowed 
points lie on the 

~ /i = 1 TeV (9) 

contour. For the chosen value of the top mass, rrit = 178 GeV, the DM constraint (jHI) is 
satisfied on this contour for mi/2 > 3 TeV and mo > 6.5 TeV; for the new preliminary world 
average top mass of 173 GeV 1^^, the lower bound on mo would be reduced to ~ 5.5 TeV. 
This is the higgsino LSP domain of mSUGRA (gaugino fraction Zg ^ 0.1).^ In this region 
there is a near degeneracy among the lighter chargino and neutralino states, i.e. 

m^ ~ m^o ~ m-+ ^ /i ^ 1 TeV. (10) 
So the major annihilation processes correspond to the pair and co-annihilation 

wwizz), m ^ hh, ^ //, (11) 

where i = 1,2 and j = 2 (1) for i = 1 (2). Although Z couples to a pair of x via their 
higgsino components, the coupling is proportional to the difference N^^ — Nf^ jTTj . Hence it 
vanishes in the limit of Mi, |/i| ^ Mz, where the Xi 2 eigenstates correspond to the symmetric 
and antisymmetric combinations of and ifg- In the same limit, the off-diagonal ZX1X2 
coupling reaches its maximal value. Thanks to the annihilation processes fllll) . the string of 
points satisfying the constraint (jS)) continues indefinitely upwards on the = 1 TeV contour, 

^There is also a small allowed region ,12 with ~ rnh/2, h being the lighter CP-even Higgs boson; 
this is however not visible at the scales chosen in Fig. 1. 

^There must be allowed points also in between the focus point and TeV higgsino-LSP regions. However, 
this DM-allowed strip is very narrow, since the transition between a lighter higgsino as LSP, with too small 
a relic density, and Bino-like LSP with much too high a relic density, is very rapid. The scan of parameter 
space used in Figs. 1 therefore found no allowed points for toq between 4 and 6.5 TeV. 
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whereas all other DM-allowed regions in mSUGRA are finite in the (^«i/2, '"^o) plane. In this 
sense the constraint (jHI) favors the higgsino LSP domain of the mSUGRA model. 

On the other hand, Figs. 1 imply that in this domain all superparticles are quite heavy. 
We just saw that even the LSP has a mass near 1 TeV. Moreover, we needed mi/2 > 3 TeV 
and tuq > 6.5 TeV for rrit = 178 GeV. This means that the electroweak gaugino {B, W) 
masses are at least in the few TeV range, while the masses of gluinos and all scalars (except 
for the lightest Higgs boson) are near 10 TeV or even higher. This aggravates the "little 
hierarchy" problem considerably; however, the finetuning required is still very much smaller 
than in split supersymmetry [A^. 

Higgsino and gaugino masses in the above range are still compatible with gauge coupling 
unification within the uncertainty of GUT scale thresholds. A sfermion mass scale near 10 
TeV is adequate to solve the problems of fiavor and CP violation even without assuming 
fiavor universality fHl- This leads us to more general SUSY models, which we will comment 
on in Sec. 5. In the next two sections we investigate the prospects of probing scenarios 
with heavy higgsino-like LSP in collider and dark matter experiments. These phenomeno- 
logical investigations are largely model-independent, so long as the remaining sparticles lie 
significantly above the higgsino-like states. 

3. Probing the Higgsino LSP Region in Collider Exper- 
iments 

Sfermion and gluino masses of ^ 10 TeV and electroweak gaugino masses of at least a few 
TeV put them well out of reach of the LHC. The only superparticles which can be produced 
there with significant rates are the nearly degenerate charged and neutral higgsinos, xt ^i-nd 
Xi 2 of mass ~ 1 TeV. We have computed the mass differences including radiative corrections 
[T7] and found them to be restricted to the range 

Srric = — m^o < 10 GeV, (12) 

with m^o — m^o ~ 2Smc for tan/5 ^ 1. Thus the xf and X2 decay products will be too soft 
to be detected efficiently on top of the underlying event at a hadron collider. Therefore one 
has to tag the pair production of higgsinos at the LHC. The by far best tag is provided by 
the two forward jets j in x pah production via vector boson fusion (19ij . 

PP xtx^jj, xtxiJJ, XiX2jj = 1, 2). (13) 

We have computed the resulting higgsino signal at the LHC closely following [19j and a similar 
investigation for an invisibly decaying Higgs signal in [201 ■ The selection criteria used are: (i) 
two forward jets in opposite hemispheres with Ej, > 40 GeV and 2 < \rij\ < 5; (ii) Arjjj > 4; 
(iii) Mi^^{jj) > 1200 GeV; (iv) > 100 GeV; (v) Acpjj < 57° and (vi) the central jet veto 
as defined in ref. ||2^. The backgrounds come from vv) and W{-^ iu) production 

via electroweak (vector boson fusion) and QCD (higher order Drell-Yan) processes where 
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C is assumed to escape detection for < 10 GeV. Following PD] we have assumed the 
efficiency of the central jet veto to be 0.9, 0.82 and 0.28 for the signal, electroweak and 
QCD backgrounds respectively. The total background is 64 fb assuming conservatively the 
renormalization scale for as to be the lower jet--Er. One expects to measure the background 
to a high precision of ~ 1.2% from the visible Z ^ II and W —>■ iu events [201 • Adding 
this uncertainty in quadrature to the statistical error on the background, and assuming an 
integrated luminosity of 100 fb~^, one thus needs a signal cross section of at least 5.5 fb for 
a 5(7 discovery. Unfortunately the cross section after cuts for the production of higgsino-like 
charginos and neutralinos with mass near 1 TeV is several orders of magnitude below this 
value. We therefore conclude that the LHC will not be able to probe the region of parameter 
space we are interested in. 

The most promising machine for detecting a 1 TeV higgsino LSP is the proposed 3 TeV 
linear e+e~ collider CLIC j2I] . We shall follow the strategy of ref. [S^ for computing the signal 
and background cross sections. The same strategy has been followed by the LEP experiments 
for setting mass limits on a higgsino LSP in particular the OPAL collaboration has 
used it to set a mass limit of 90 GeV. The higgsino pair production is tagged by a photon 
from initial-state radiation (ISR), i.e.^ 

e+e" ixtxl, 7X1X2 ■ (14) 

If the xf and X2 decay products remain undetected, the main physics background is 

e~^e~ -^uu. (15) 

The photon is required to have an angle 9 > 10° relative the beam axis. Moreover it is 
required to satisfy 

7717 ^ 7717 min /— Sin^min /Tr\ 

1 + Sm 6'min 

which vetos the radiative Bhabha background e+e^ 7e^e^, by kinematically forcing one 
of the energetic to emerge at an angle > 6'min- At CLIC energy of = 3 TeV, 

"^"^ = 50 (100) GeV for ^.^i^ = 1°(2°). (17) 

The OPAL detector has instrumentation for detection down to ^^min = 2°, while it seems 
feasible to have it down to 1° at the future linear colliders [22] • We shall show results for 
both Ej^ ™™ = 50 and 100 GeV. We shall also impose the recoil mass cut 

Mrec = v^(l - 2E^/v^)i/2 > 2m^, (18) 

which is automatically satisfied by the signal ()14j). Fake photon background processes have 
been effectively suppressed by the OPAL collaboration [22] by requiring photon isolation and 
a minimum value for the total px-, which are automatically satisfied by the signal as well as 
the e^e^ '-fuu background. Therefore we shall not impose these requirements. 

■^The cross sections for Xi = 1)2) pair production are negligible, since the diagonal XiXi^ couplings 
are very small, as remarked earlier. 
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Figure 2: Cross sections for the higgsino signal ()14|) and the neutrino background (jlSj) at 
CLIC {y/s = 3) TeV, produced with a photon tag of Ej, > 50 GeV. Initial state radiation is 
included. 

Fig. 2 shows the signal and the background cross sections from (fT^ and ()15|) respectively 
against the higgsino LSP mass, for = 50 GeV. In calculating these cross sections, we 

have included initial state radiation (ISR) effects by convoluting the hard 2^3 cross sections 
with electron distribution functions, as described in ref.|21]. This allows background events 
with on-shell Z boson, if there is an energetic ISR photon going down one of the beam pipes; 
ISR therefore increases the total background by a few %. We also computed the higher-order 
background process 

e^e^ Zz/z/7, where Z uu, (19) 

and found it to contribute about 10 fb after cuts - much less than the background (|TH|) . but 
still significantly more than the signal. 

The signal cross section is reduced by ISR by ~ 10% for = 1 TeV, since it effec- 
tively reduces the amount of phase space available. The same effect increases the signal for 
smaller LSP mass, since it increases the s— channel photon and Z propagators. The signal 
is dominated by the chargino pair production. For 1 TeV higgsino mass one expects a signal 
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cross section of only ~ 0.8 fb against a background of ~ 1050 fb. Thus for the projected 
CLIC luminosity of 1000 fb~^ one expects 800 signal events against a background of 10®, 
corresponding to Ns/y^Ns — 0.8 only. Evidently it is a hopeless situation unless one can 
suppress the background (fT3j) by identifying the soft xf and X2 decay products. This remains 
true when the cut on Ej^ is increased to 100 GeV (not shown). 

The method of identifying these particles would depend on the decay length cr, which 
depends strongly on the mass difference 6171^. This decay length has been estimated in ref. 
[22] ioT a specific model of an iso-triplet chargino. It is shown there that for 6mc < 1 GeV 
one expects to detect the chargino track and/or a decay 7r^(£^) track with displaced vertex in 
a standard micro-vertex detector. One can easily check that the decay length of the charged 
higgsino is about twice as large as the iso-triplet chargino of [22] • Hence these tracks should 
be even more clearly detectable in this case. 

But one expects prompt chargino decay for Srric > 1 GeV, which holds over most of our 
parameter space of interest. For this case the OPAL collaboration j23| has found that the 
resulting charged tracks can be detected with > 50% efficiency for the signal, and used it 
to eliminate the juu background. For the present case such an efficiency corresponds to 
a respectable signal size of ^ 400 events. However, a new problem arises at future linear 
colliders, which did not occur at LEP. The large charge density in the bunches gives rise to 
"beamstrahlung" when the two bunches cross. The collision of beamstrahlung photons can 
then form an underlying event containing several soft particles j2S]. If this happens in the 
same bunch crossing as a hard e~^e~ '-/uu annihilation, one obtains a similar final state 
topology as in the signal. 

It is not possible to speculate at this stage on the level of this underlying event back- 
ground at CLIC. All we can say is that an underlying event resembling the xt xf decay 
products must not occur in more than 1% of all bunch crossings. Neglecting efficiencies, this 
would correspond to ~ 10^ background against ~ 800 signal events for = 1 TeV - i.e. 
Ns/\/Nb ~ 8. We will see below that one can tolerate a higher level of underlying events 
from beamstrahlung if the e~^e~ beams are polarized. Finally, we note that beamstrahlung 
will also change the effective beam spectra, and hence the cross sections ()14j) and p5|l . 
These effects will need to be included once the beam characteristics have been fixed. 

The reason for the large cross section for the background ()15p compared to the signal 
processes (fT^ is the t— channel W exchange contribution to the background. This can be 
suppressed with right- (left-)handed polarization of the e~ (e"*") beam. In fact it is easy 
to see that for 100% polarization of one of the beams the background cross section will go 
down to the level of the signal. This is not feasible, of course. What we shall do instead is 
to estimate the signal and background for the same beam polarizations as envisaged for the 
ILC i-e- 

Pg- = 0.8 (mostly right — handed) and Pe+ = —0.6 (mostly left — handed). (20) 
It is easy to check that this corresponds to the following fractional luminosities 

enet : ele+ : eZe+ : e^e+ = 0.72 : 0.02 : 0.08 : 0.18, (21) 
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while each was 0.25 in the unpolarized case. The dominant contribution to the background 
(|15|) from t-channel W exchange contributes only to the second combination e^e^. Hence it 
is suppressed by a factor of .02/. 25 = .08. The higher-order background (fT^ is suppressed 
by a similar factor. One can also check that the xtxi and X1X2 contributions to the signal 
(fT^ are modified by factors of 0.6 and 1.3 respectively, resulting in an overall suppression of 
the total signal by a factor 0.8. Fig. 3 shows the total signal and background cross sections 
for Ej, > 50 and 100 GeV. In either case one gets a Ns / \J Nb — 2.1 for the CLIC luminosity 
of 1000 fb-^ But one has a better Ns/Nb ~ 400/39, 000 events for the E:^ > 100 GeV cut. 
Recall that this cut requires instrumentation down to 2° instead of 1° to eliminate the 7e"^e~ 
background. Hence this harder cut seems advantageous to us. 
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Figure 3: Cross sections of the higgsino signal ()14|1 and neutrino background ()15|1 at CLIC 
with polarized e~ and beams. Initial state radiation is included. 

Fig. 4 shows the recoil mass distribution of the background (fT^ along with that of 
a 1 TeV higgsino signal ()14|) for both Ej. cuts and polarized beams; ISR has again been 
included. As discussed above the background can be suppressed if the soft xf or X2 decay 
products can be detected. In case of the background, such soft particles can only come from 
beamstrahlung reactions like 77 — > 7v~^7r~ ,£~^i~ , ■ ■ ■ underlying the background (fTH|) . Not all 
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such reactions will lead to events with similar characteristics as the signal. For example, 
one can envision applying cuts on the angular distribution of the soft particles, which tend 
to peak at small angles in two-photon events, but are quite central for most signal events. 
Another possible discriminator is the pt imbalance of the soft particles (i.e., not counting the 
hard tagging photon), which is expected to be larger for the signal than for the background. 
In devising such cuts, the characteristics of the (largely non-perturbative) background can 
be taken from measurements in the pure background region M^-^c < 2m^. Comparing the 
observed cross section for the remaining background over the Mrec < 2 TeV region with 
the prediction of Fig. 4 would give an estimate of the fraction of surviving background due 
to beamstrahlung. Since this fraction should be independent of Mrec, one can use this to 
estimate the cross section of the surviving background in the Mfgc > 2 TeV signal region. 
Any excess over this estimate would represent the higgsino signal. One might even be able 
to estimate the higgsino mass from the threshold of the excess cross section. It is easy to 
see that underlying events from beamstrahlung at the 10% level (after cuts) correspond to 
a reduction of the background to 10% and hence will increase Ns / y/ Nb ratio from 2 to ~ 6. 
Thus with polarized beams one can tolerate the underlying event at the 10% level. 
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Figure 4: The recoil mass distributions of a 1 TeV higgsino signal ()14j) and the neutrino 
background (fTH|l at CLIC with polarized e~ and e"*" beams. Initial state radiation is included. 
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4. Higgsino LSP Search in DM experiments 



One can see from the second paper of ref. |HI that the higgsino LSP signal is too small 
to be measurable in direct dark matter search experiments. The reason is that the signal 
comes from spin-independent XP scattering, which is dominated by the Higgs boson {h, H) 
exchange. Since its coupling to a x pair is proportional to the product of their higgsino and 
gaugino components, it is very small for a higgsino dominated LSP. The signal is further 
suppressed by the large LSP mass. 

We have also checked that the neutrino signal coming from the x P^i^' annihilation at the 
solar core is too small to be measurable at an IceCUBE size detector. Here the signal size is 
determined by the spin-dependent xP scattering cross section via Z boson exchange, which 
is very small due to the suppressed diagonal Zxx coupling; see the remark following eq. (fTT|) . 

The most promising signal for TeV higgsino DM comes from the pair annihilation pro- 
cesses 

XX 77, XX iZ, (22) 



resulting in a monochromatic 7— ray line p7t 128]. The dominant contributions to these pro- 
cesses come from W^Xi loops, and are suppressed by only a factor in the denominator 
instead of m^. This results in a large cross section for 7— ray production from (j^^ for a TeV 
scale higgsino, 

va^^ ~ va^z ~ 10"^^ cm^s"\ (23) 

where v is the velocity of the DM particles in their cms frame. The resulting gamma ray 
flux"^ coming from an angle ip relative to the galactic center is given by 

A, 

where p{i) is the dark matter energy density and = 2 (1) for the 77 {'jZ) production 
process. This can be rewritten as fI7\ 



N^va \ (I TeV\ 



2 



m) = 1.87 X 10-" (^jrst^j [— ) J(*)cm-s-'s.-' , (25) 
where 

J{^) = [ p^{t)M{'4j)/ [(0.3 GeV/cm=^)2 ■ 8.5 kpcj (26) 

J line of sight '- 

is the line integral scaled by the squared DM mass density in our neighborhood and by our 
distance from the galactic center. 

Several Atmospheric Cerenkov Telescopes (ACT) have started recording or are on their 
way to record such 7— rays from the galactic center - i.e. MAGIC and VERITAS in the 



^It has been pointed out very recently |20j that tree-level higher order processes, in particular XX ~* 
W~^W~j, can increase the flux of photons with ~ by up to a factor of 2. While significant, this 
enhancement is still much smaller than the uncertainty coming from the DM distribution near the center of 
the galaxy, as discussed below. 
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northern hemisphere and H.E.S.S. and CANGAROO in the south. One generally expects 
a concentration of DM in the galactic center; but its magnitude has a large uncertainty 
depending on the assumed profile of the DM halo density distribution [HDl EH E2]- The 
cuspy NFW profile corresponds to 

(J(0))ac=.ooi ^ 1000, (27) 



which represents the DM flux in the direction of the galactic center averaged over the typical 
ACT aperture of AQ = .001 sr. Extreme distributions, like the spiked profile [HI] and core 
profile correspond to increase and decrease of this flux respectively by a factor of ~ 10^. 




Figure 5: Monochromatic gamma ray flux from DM pair annihilation near the galactic center 
shown for the NFW profile of DM halo distribution and an aperture AQ = 10"^ sr. 

We have computed the 7— ray line signal (j25p for the NFW profile and the aperture 
AQ = .001 sr using the Dark SUSY code |33 . Fig. 5 shows the resulting signal against the 
DM mass, where we have added the 77 and 7Z contributions, since they give identical photon 
energy (= m^) within the experimental resolution. This result agrees well with that of ref. 
j34j . The vertical spread in the higgsino band reflects the dependence of the annihilation 
cross section (j23|) on the mass difference 6mc. As noted in |34j the discovery limit of the 
above-mentioned ACT experiments goes down to 10~^^ cm~^s~^. Thus for the NFW profile 
one expects a 7— ray line signal that should be detectable for the WMAP favored mass range 
of ~ 1 TeV. Recall that the signal rate will go up by a factor or of ~ 10^ for the spiked 
profile [SI], while it will go down by a similar factor for the core profile [S2]- In the latter 
case it will fall below the discovery limit of these ACT experiments. 

In fact, H.E.S.S. did detect TeV photons coming from the direction of the galactic center 
|35j . However, they observe a continuous spectrum extending beyond 10 TeV in energy, 
which can be described quite well by a power law. This is not what one expects from DM 
annihilation. In fact, the spectrum looks very similar to that of other "cosmic accelerators" 
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observed by the H.E.S.S. telescopes. It is currently not clear whether this signals comes right 
from the center of our galaxy (defined as the location of a supermassive black hole), or from 
a nearby supernova remnant (SNR); in particular the SNR Sagittarius A East might be the 
culprit [36j. Note that SNR are known to emit TeV photons with power-law spectra. The 
H.E.S.S. collaboration is now working on improving their angular resolution. If the source 
of the observed TeV photons, whose flux is well above the detectable limit, continues to 
coincide with the galactic center within the resolution, the discovery of a line signal from 
DM annihilation at the center would become more difficult, since one would then have to 
look for a peak in the spectrum on top of a sizable smooth background. 

5. Higgsino LSP in other SUSY Models 

We have so far concentrated in the mSUGRA model for its simplicity and economy of 
parameters. However one can easily see that all our results hold for the higgsino LSP in 
a host of other SUSY models. This is because the relevant interactions are the higgsino 
interaction with the gauge bosons, which are completely determined by the gauge charges 
of Hi and H2 along with their mixing. Both these features are common to all variants 
of the minimal supersymmetric standard model (MSSM). Thus the cross sections for the 
higgsino annihilation processes (jlip . and the resulting higgsino masses (fTUI) obtained using 
the constraint (jH)) on their relic density, are common to all MSSMs, as long as the other 
superparticles and heavy Higgs bosons have masses ^ 2 TeV. The same is true for the 
signals depicted in Figs. 2-5. This also explains why the DM results of Figs. 1 and 5 
are essentially independent of all SUSY parameters except for the higgsino masses. We 
therefore have made no mention of these parameters while presenting the collider signals 
of Figs. 2-4. It should be mentioned here that the strip to the left of the DM allowed 
higgsino LSP range in Fig. 1 corresponds to an underabundance of DM relic density in the 
standard cosmological model. However, additional thermal and nonthermal mechanisms of 
DM production have been suggested [SZlEHl, which could enhance the DM relic density over 
its standard cosmological model value. In the presence of such mechanisms the DM allowed 
range will move to the n < 1 TeV region; and so will the higgsino LSP mass. In that case 
the collider and DM signals shown in Figs. 2,3 and 5 over the LSP mass range of 200 - 1000 
GeV will become relevant. 

We saw above that in the context of mSUGRA, a TeV higgsino can be the LSP only 
if sfermions lie near 10 TeV or even higher. This is adequate for suppressing FCNC pro- 
cesses even without assuming flavor universality of scalar soft breaking parameters. It also 
allows 0{1) phases in the soft breaking sector without violating constraints on CP violating 
processes. This greatly opens up the allowed parameter space, even if one keeps the two 
soft breaking Higgs masses the same in order to achieve radiative EWSB. Examples for such 
models are the so-called inverted hierarchy and more minimal supersymmetry models |18] . 
At the cost of additional finetuning(s) one can even entertain the idea of moving the 
sfermion masses to yet larger values, as in the split SUSY model [3]. 

We saw in Figs. 1 that mSUGRA predicts the LSP to be Bino-like over most of the 
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theoretically allowed parameter space. This can be traced back to the fact that the coefficient 
C2 in eq.© is quite large and positive, while |Ci| is small, so that > Mi at the weak 
scale unless ttIq ^ '^\/2- C*i can be increased if the Higgs soft breaking masses exceed the 
stop masses at the GUT scale [101 • On the other hand, C2 can be reduced if the ratios 
M1/M3 and/or M2/M3 are increased IH] relative to their mSUGRA values Q. Models with 
non-universal scalar and/or non-universal gaugino masses therefore often can accommodate 
a higgsino-like LSP more easily than mSUGRA does. Finally, if one reduces the input 
scale, i.e. the scale where supersymmetry breaking is mediated to the visible sector one 
simultaneously increases Ci and reduces C2, again making it easier to obtain a higgsino-like 
LSP. All our results will apply equally to these models. 

6. Summary and Conclusions 

We have seen that a higgsino-like LSP can be Dark Matter in a variety of supersymmetric 
models. In the most constrained case, the mSUGRA model, this remains a possibility for 
arbitrarily large values of mo and mi/2, thereby greatly enlarging the cosmo logically allowed 
region of parameter space. As discussed in Sec. 5. a higgsino-like LSP can also be realized 
in many extensions of the mSUGRA model. 

In standard cosmology, and assuming that the LSP was in thermal equilibrium after the 
period of last entropy production, the LSP relic density can be calculated uniquely from its 
( CO-) annihilation cross sections. In case of higgsino-like LSP one finds that a mass near 1 
TeV is required. This makes sparticle searches at colliders quite challenging. In most models 
strongly interacting sparticles have masses at least a factor of 5 above the LSP mass; this 
is true in particular for all models with (approximate) gaugino mass unification near the 
scale of Grand Unification. This means that the usual SUSY signatures at the LHC will not 
work. We found in Sec. 3 that the production of two higgsino-like states in vector boson 
fusion also does not give rise to a detectable signal at the LHC if these states lie near 1 TeV. 
Moreover, the energy of the next (international) linear collider ILC will not be sufficient to 
produce pairs of TeV sparticles. 

We therefore have to consider more futuristic colliders. We saw in Sec. 3 that the pro- 
posed 3 TeV e^e~ collider CLIC offers quite good prospects, if the level of beamstrahlung 
induced underlying events can be kept under control. This can be achieved by designing 
the accelerator such that the flux of beamstrahlung photons remains small, and/or by build- 
ing a sufficiently sophisticated detector so that the kinematic distributions of soft particles 
produced in two-photons events can be distinguished from those of the soft decay products 
of the heavier higgsino-like states X2 ^^^1 xf • We also saw that the ability to polarize the 
incident beams would be very helpful. 

In order to show that a given particle forms the Dark Matter in the universe, one will 
eventually have to detect these relics. We saw in Sec. 4 that in case of a higgsino-like LSP 
the most promising search is that for a 7— ray line at E^i ^ m^. The fiux of such photons 
should peak in directions where DM particles accumulate. The by far most promising site 
is therefore the center of our galaxy. Unfortunately here the signal might be masked by 
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the recently observed flux of TeV photons with a continuous spectrum extending beyond 10 
TeV. Improved angular and/or energy resolution would be helpful in enhancing the signal 
to background ratio in this case. 

We conclude that a TeV higgsino is a viable supersymmetric Dark Matter candidate. The 
large sparticle masses characteristic of such a scenario require some amount of finetuning, 
but alleviate problems with flavor-changing neutral currents and CP violation. Testing 
this scenario experimentally is challenging, but should be possible at future multi-TeV e^e~ 
colliders like CLIC, and perhaps through the observation of a TeV 7— ray line in atmospheric 
Cerenkov telescopes. Finding TeV higgsinos either at colliders or in Dark Matter search 
experiments is certainly easier than finding gravitinos, which have been much discussed 
lately as possible Dark Matter candidates This scenario should therefore be taken 

seriously, in particular if the LHC fails to discover super symmetry. 
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